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Abstract 

We classify the gapped phases of Zjv parafermions in one dimension and construct a represen- 
tative of each phase. Even in the absence of additional symmetries besides parafermionic parity, 
parafermions may be realized in a variety of phases, one for each divisor n of N. The phases can 
be characterized by spontaneous symmetry breaking, topology, or a mixture of the two. Purely 
topological phases arise if n is a unitary divisor, i.e. if n and N/n are co-prime. Our analysis is 
based on the explicit realization of all symmetry broken gapped phases in the dual Z^v-invariant 
quantum spin chains. 

1 Introduction 

Topological phases of matter received a lot of interest recently due to their peculiar properties such as 
the existence of robust edge modes, non-abelian excitations and their potential applications to quantum 
computation [HE]- While one-dimensional bosonic systems do not exhibit intrinsic topological order 
[3], the situation is very different for fermions. According to the general theory, there should exist two 
phases of fermions which are described by a Z2 topological invariant [HEJE]. An explicit realization of 
both phases is provided by Kitaev's Majorana chain which models the dynamics of spinless fermions 
in a quantum wire, in proximity with a p-wave superconductor [7J. 

The most remarkable feature of Kitaev's model is the emergence of stable gapless Majorana edge 
modes in the topologically non-trivial phase [7J. In view of their non-abelian statistics, being able 
to create and to manipulate such isolated Majorana fermions would be a major step towards the 
realization of topological quantum computers [SJ. While Kitaev's chain merely serves as a proof 
of principle, by now more realistic experimental setups for the creation and detection of Majorana 
fermions have been suggested [TTJ] . The basic ingredients are spin orbit quantum wires in proximity 
with an s-wave superconductor and in a magnetic field. Recent experiments based on this setup 
provide first evidence that Majorana fermions indeed exist |11[ [T2] even though at least part of the 
observed features could also be explained by disorder effects [T3J HH [15] . 

More recently, it has been investigated whether similar phenomena can be realized in heterostruc- 
tures involving fractional topological insulators. It was found that such systems can be used to isolate 
parafermions which may be regarded as fractionalized Majorana fermions |16[ 117] , These efforts re- 
ceived their motivation from notes of P. Fendley on the existence of edge modes in parafermionic 
chains, unpublished for a long time but eventually made available in Ref. |18j . Other work on this 
subject includes Refs. [T9l [20] . 

The study of parafermions has a long history. The exploration of para-particles began in Refs. |21|, 
[22] in an attempt to generalize the canonical commutation relations of bosons and fermions. In 
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contrast, these days the predominant notion of parafermions arises in the context of statistical physics, 
in connection with the so-called Z at clock model [23J . This model generalizes the Z2 symmetry of the 
Ising model in a natural way to the Z^ symmetry of a clock with N marksQ 

Just as the Ising model may be mapped to a problem of free fermions, using a variation of the 
non-local Jordan- Wigner transformation, the Z^r clock model exhibits a dual description in terms 
of fundamental parafermionic degrees of freedom. The parafermionic operators satisfy a Zj\r-graded 
variant of the usual Clifford algebra relations. This duality generalizes the well-known duality between 
the Ising model and free fermions for N = 2. It is known that the topologically trivial and the non- 
trivial phase of the Majorana chain are equivalent to the ordered and the disordered phase of the Ising 
model, respectively j26j. 

The previous paragraphs suggest two natural avenues to the classification of parafermionic chains. 
One first of all can generalize the group theoretical approach of Ref. [6] that has been applied suc- 
cessfully to enumerate and characterize the different phases of Majorana chains with or without time- 
reversal symmetry. This approach has the advantage of yielding the maximal set of possibilities in a 
clean and conceptually clear way. On the other hand, the proof that all topological phases can actually 
be realized requires a method of combining phases such that their topological charges add up. This 
procedure might not necessarily lead to the most natural representative of the new phase. Moreover, 
the approach is limited to phases without broken symmetry. This approach has been followed recently 
in Ref. 1361. 



The second approach which shall be followed in this paper is a complete characterization of gapped 
symmetric and symmetry broken phases of Z^r quantum spin chains. These can be engineered system- 
atically on the level of ground states and their associated frustration free parent Hamiltonians |27| . 
Via the duality mapping mentioned before, this classification provides an explicit characterization of 
all parafermionic phases which cannot be connected without crossing a phase transition, i.e. closing 
the gap. 

We shall find that the distinct parafermionic phases are labeled by the possible divisors of N. 
Generically they exhibit coexistence of both ferromagnetic and topological order. Only part of the 
degeneracy of the ground state of an open chain is topologically protected and due to gapless edge 
modes, and the rest can be lifted by spontaneous symmetry breaking. When the chain is closed, 
only the topological degeneracy disappears, leaving in general with a degenerate ground state. In 
particular, for a given divisor n of N, the phase will be purely topological if N and N/n are co-prime, 
and will exhibit a unique ground state in a closed system and gapless edge modes in an open system. If 
instead N/n = gcd(n, N/n) the phase is topologically trivial and entirely characterized by symmetry 
breaking. The details of the classification will be explained through concrete realizations of the phases 
in the text. 

We would like to stress that our classification distinguishes only phases which cannot be connected 
by any local Z^r-invariant interaction. Clearly, specific models with an explicit Hamiltonian depending 
on a given set of coupling constants may have arbitrarily complicated phase diagrams. 

The paper is organized as follows. In Section [2] we introduce the kinematical setup that we are 
working with. This includes a detailed discussion of the symmetry algebra and of the Hilbert spaces 
in question. The possible symmetry breaking patterns and order parameters of the Z^r quantum 
spin chain are analyzed in Section [3] and a representative Hamiltonian is constructed for each phase. 
These results form the heart of our paper. Finally, we perform a non-local duality transformation in 
Section [3] and interpret our previous results from a parafermionic perspective. As it turns out, the 
original symmetry broken phases of the spin chains generally exhibit a mixture of symmetry breaking 

x The name "parafermions" is sometimes also used for the basic fields in a family of Zjv symmetric conformal field 
theories that arise at critical points of the clock models [24]. These parafermions may be used to construct the Read- 
Rezayi states, i.e. they are of relevance for the study of fractional quantum Hall systems [25]. In the present paper, we 
shall be concerned exclusively with the gapped phases of Zm parafermions. 
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and topological protection. We end with concluding remarks and an outlook to open problems. 



2 Parafermions and Z^v spin chains 

We introduce the mathematical framework within which our considerations take place. This includes 
a detailed discussion of the symmetry algebra and of the Hilbert spaces involved. 



2.1 



spin chains 



Zjv quantum spin chains are Zjy symmetric models where the spins can show in any of N directions 
in the plane which can be thought of as corresponding to the iV distinct N th roots of unity on the 
unit circle. This setup generalizes the familiar Z2 case of the Ising model where one deals with spins 
pointing either up or down, and at the same time, provides a specialization of the XY model in which 
the spins are restricted to take discrete values. 

If we denote the size of the chain by L, the Hilbert space of the chain will be given by 

(C N )® L . 



We choose a basis in H given by \q±, . 
operators 
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and as obvious from the definition, di,Tj commute if % 7^ j. Note that for N = 2: a = a 2 , r = a x . 
The space of linear operators on % is spanned by monomials Yli=i a T T T ' ■> Pi'Qi ^ ^n- 

The global symmetry G = Z^r is generated by the operator shifting the labels by one unit on every 



site: 



P 



n 



(2.7) 



and each tensorand C N in the Hilbert space T~L is direct sum of iV one-dimensional irreducible repre- 
sentations of Z^v- The Hamiltonians H we are going to consider will be invariant under the action of 
Zjy. In other words, [H, P] = 0. 



2 These are the relations of the Heisenberg group which is defined by exponentiation of the Weyl algebra [Q,P] = 
QP — PQ = il. In terms of the generators a — exp(iP), r = exp(iQ) and the central element Z = exp(il) the group 
multiplication reduces to <j t t s = Z rs T s a r for any parameters r,s£l. 
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2.2 Jordan Wigner transformation 



Now we perform a generalized Jordan Wigner transformation to define the parafermions [23]. First 
we define operators living on the dual lattice as: 



Mi+l/2 = > ft+l/2 = II r fc ' 



fc=i 

so that Tj = ft-i/2Aj + i/2' an< ^ Z 2 '' 3 sa fisfy the same commutation relations as respectively r, a. p is 
the disorder operator dual to the order operator a, implementing the action of symmetries on spins 
between sites 1 and i, namely shifting those spins by one unit. Parafermions are then defined on 
middle points between vertices of the original and dual lattice, as the product of order and disorder 
variables: 

72i-i = (j[ Oi , 72i = J N ^)I 2 Qj r^j a in . (2.9) 

They satisfy the following commutation relations 

7 f = l, 7i7i = <«>7j7< (for l<i<j<2L), (2.10) 
and the definition can be inverted to yield: 

wU-uP+W-Wi&l (2-11) 

n = u-^l^_ ll2l . (2.12) 

7's are the generators of an algebra which generalize the Clifford algebra of "Majorana fermions" 
occurring at N = 2 [28]. Elements of a generalized Clifford algebra (with even number of generators, 
as in the case at hand, where the algebra is isomorphic to the complex matrices acting on the Hilbert 
space of the chain) are the sum of homogeneous elements X which have a well defined Z^-grading 
|X| G Zjv given by: 

FXP^ = J%, (2.13) 
with P the symmetry generator expressed in terms of the parafermions as 

P = n^- 1)/a ^7iii. (2-14) 
i=i 

In particular, the parafermionic generators themselves are elements of degree 1. 



3 Symmetry broken phases in Z^v spin chains 

The general theory developed in [271 [29] classifies gapped phases of spin chains with symmetry group G 
in terms of classes of projective representations of subgroups H to which the symmetry is spontaneously 
broken. The subgroups of Z^v are given by Z n , where n is a divisor of N. Z^r only has trivial projective 
representations for any TV (see [30l Appendix J]), so that gapped phases of Z^r spin chains are indexed 
by divisors of TV and fully characterized by spontaneous symmetry breaking (SSB). 

Statistical systems with Zj\r symmetry have been studied extensively over the years, and their 
phase diagrams have been discussed in a variety of places, see e.g. [311 |32l [231 [33] ■ 111 this section, 
we explicitly construct representatives of all gapped phases we expect to have. On the way, we shall 
emphasize important aspects related to order parameters which will play a key role in the description 
of parafermionic phases in Section [H 
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3.1 Construction of the phases 



We first briefly recall the defining features of SSB in Z^r spin chains. The breaking of Z2 symmetry 
in the ordered phase of the quantum Ising chain is an elementary textbook example |34j : the phase is 
characterized by a non-vanishing expectation value of the order parameter a. The generalization to 
the breaking of Z^r down to Z n , n divisor of N, is straightforward. 

In the thermodynamic limit, the local order parameters describing the properties of the spin chain 
are given by the powers a an with a € Z m — {0} and m = N/n. The residual symmetry which leaves the 
order parameter invariant is Z n . It is generated by P m , the m th power of the global parity. (If n = N 
the phase is disordered.) Elements of Zjy/Z n = Z m are generated by P n , and permute the ground 
states cyclically among themselves, so that the ground state manifold is m-fold degenerate. We will 
refer to Z m as the broken symmetries. The degeneracy of the ground state encodes a permutation 
representation of the symmetries on the ground state manifold. This representation corresponds to 
an invariant of the phase and cannot change unless a phase transition occurs |27J. The symmetry 
breaking perturbation h^2Aa? n + h.c.) (h small and positive) of the Hamiltonian will select a single 
ordered ground state, leading to a non-zero expectation value of the order parameter even when the 
external field h is set to zero. 

For later convenience we will now explicitly construct representatives of each phase. This includes 
both the desired states and an associated parent Hamiltonian. According to the above discussion, a 
phase labeled by n should have m ordered ground states \ip a ) (a G Z m ) which satisfy 

Ptya) = \^ a+1 ) , P km \^ a ) = \il> a ) , (3.1) 

for all k £ Z n . A simple solution of these constraints is given by the product states 

L n-1 

\ipa) = (^) , where \a)) : = \a + mk) . (3.2) 

i=l k=0 

The states \a)) are orbits under the subgroup Z n of Z^r. As such they span the subspace C m C C N 
and they satisfy the periodicity condition \a)) = \a + m)). 

The associated Hamiltonian will be defined in such a way that the m-fold degeneracy of the ground 
state is already exact in finite systems and not only in the thermodynamic limit. This is achieved by 
choosing a frustration free combination of specific projectors which are localized on two neighboring 
sites. Taking H = — the following choice of the two-body term /ij^+i guarantees the Z^v- 

invariance and the projector property (hj i+l = /ij^+i, h\ i+1 = hi t i+i) 

- m— 1 

hi,i+i = — ^ \a))i\a))i + i((a\i((a\ i+1 . (3.3) 

The Hamiltonian H is a sum of commuting projectors, so that it is possible to diagonalize all terms 
simultaneously. The resulting spectrum is given by the integers between — L and 0, and hence the 
chain is gapped as desired. Assuming periodic or open boundary conditions will make no difference for 
the conclusions we want to draw in this section. In both cases the space of ground states is spanned 
by the vectors \ip a ) with a G Z m , for which h^i+i = 1. 

It remains to find a more explicit form of the Hamiltonian (|3.3|) in terms of the spin operators cr, 
and Tj. Using the single-site matrix elements 

1 m— 1 n— 1 

-(</!«» (Hp) = -<ff| £(™-*)^y>^b) , (3.4) 

0=0 j=o 
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it can easily be verified that Eq. (J3.3P reduces to 

m— 1 re— 1 

hm = £ E M+i)' 9 * E • ( 3 - 5 ) 

0=0 j,k=0 

We have thus succeeded in constructing gapped Hamiltonians realizing each possible SSB phase of Zjy 
spin chains. 

Before moving on, we would like to comment on the specific form of the Hamiltonians obtained. 
The cases n = 1 and n = N correspond, respectively, to the situations with no residual symmetry 
and full symmetry of the ground state. In these two extreme cases, the interaction simplifies to the 
expressions 

1 

(n, m) = (1, N) : h i}i+1 = - E M-flf (3-6) 

Q = 
1 JV_1 

(n, m) = (N, 1) : h i>i+l = ^ E T l T ^ " ( 3 - 7 ) 

i,*=o 

It is now evident that the general formula (|3.5p is obtained by multiplying together the Hamiltonians 
of a Z m chain whose symmetry is completely broken, 1/^X^/3=0 (^^(^l+i) 71 )^' that of a Z n chain 
without symmetry breaking, 1/ra 2 X^fci=o( r «™) J ( r i+i) fe - 

While our definition reproduces the conventional form of the ferromagnetic clock model deep in 
the symmetry broken phase |18j for n = 1, the case n = N leads to a more complicated expression 
than the one used in [18] , namely ^j^Sq 1 t£ (even though both Hamiltonians have the same symmetric 
ground state). Moreover, the duality transformations (|2.8p between the ordered and the disordered 
phase do not preserve the structure of the Hamiltonians (|3.5p . It is therefore reasonable to search 
for a slight modification of the Hamiltonian (|3.5p such that it reduces to the form of the spin chains 
considered in [18], while retaining the ground states above. This goal is achieved by considering the 
expression 

( 1 m— 1 -. re— 1 \ 

0=0 j=0 ) 

This operator is Hermitian, Z_/v invariant, and each single term is a projector commuting with all the 
other terms. Therefore the chain is still gapped after the modification, albeit with a different spectrum 
in comparison to H. 

The ground states are such that each term is 1 and H has exactly m ground states which are given 
by (|3.2p . Indeed, the terms l/^X]?=o 7 f" 1 ac * as identity on states (£)f = i {on £ Z m ). Demanding, 
in addition, that also l/ m SSTo 1 ( '* cr l+i)^ n ac * s as identity selects states with cc, = Hence H is 

another representative of the SSB phase under consideration. Note that the Hamiltonians H in (|3.8p 
admit a simple action (in the thermodynamics limit) under the duality transformations (|2,8p . The 
order-disorder duality simply swaps the phases labeled by n with those labeled by m, and at the same 
time interchanges a with p and r with /x. Nonetheless, we will stick to our original definition of the 
Hamiltonian, see Eq. (|3.5|) . in what follows. 

3.2 Order parameters 

In a macroscopic system with broken symmetry one will not observe a quantum superposition of ground 
states. Fluctuations will rather single out a unique ground state with a well-defined magnetization. In 
order to select a single ground state |^>o)> we a dd a symmetry breaking perturbation h^^a? 71 + h.c.) 
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to the Hamiltonian (|3.5p . Afterwards we take the thermodynamic limit and set h = 0. It is easy to 
check that this procedure implies a spontaneous magnetization 

(at):=(^Wtm^0, (3.9) 

and also (t^ 711 } 7^ for all k G Z n . More importantly, we also have 

fofiV^O. (3.10) 

In contrast, all other powers of a and p have vanishing expectation values. Now we expect that slight 
modifications of Hamiltonians of gapped phases cannot change abruptly the expectation values unless 
a phase transition is reached. Then Eqs. (|3.9] 13.10]) hold not only for the finely tuned Hamiltonians we 
have constructed but in a whole region of non-zero measure in the phase diagram, which constitutes 
the SSB phase. 

The mechanism of symmetry breaking is well understood for the Ising chain and it applies similarly 
to Zjy spins. The operator pi+1/2 disorders an ordered state by creating kink excitations, topological 
objects rotating all spins from site 1 to % by one unit. The disordered phase can be characterized by 
the condensation of these kinks [35J. 

We are now prepared to appreciate that a particular situation is encountered when there exist two 
integers a, k such that an = km. Indeed, under these circumstances we are able to construct another 
non-local order parameter 

(/f 1/2 ^r n - (3.11) 

(The non-zero expectation value is implied by Eqs. (|3,9] 13.10]) .) With the definition g := gcd(n, m), 
the solutions of the equation an = km can be parametrized as follows, 

a = k^ with fce{0,2 22 = 2Z fl . (3.12) 

We conclude that we have additional (non-local) order parameters whenever g > 1. Their existence 
is a direct consequence of being in a SSB phase of the Zjv spin chains. (Of course the presence of 
non-local order parameters does not mean that the phase at hand is topological, since the degeneracy 
of the ground state of the spin chain can be completely lifted through SSB.) In Section B~2] we will shed 
more light onto the physical implications of the presence of this new non-local order parameter. It is 
no coincidence that the product (|3.1ip of order and disorder variables that we have just encountered 
agrees with powers of the parafermion, see Eq. (|2.9p . 



4 Classification of parafermionic phases 

In this section we will translate the classification of symmetry breaking patterns in Zjy spin chains 
into the dual language of parafermions. It will be found that the resulting parafermionic Hamiltonians 
generally exhibit a mixture of symmetry breaking and non-trivial topology which can be made manifest 
in the ground state degeneracies and in the occurrence of protected edge modes. We shall point out 
certain subtleties that arise when considering closed boundary conditions as opposed to open ones. 

4.1 Duality transformation of the Hamiltonian 

Our first goal is to rewrite the Hamiltonian Eq. (|3.5p in terms of the parafermions. Using the trans- 
formation rules (|2.1ip and the relation 

( 7 f 7j 6 ) c = w -«*<c-l)/2^ (for {<j ^ 
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we immediately obtain the two-site Hamiltonian 

m— 1 n— 1 

/3n(iV+/3n)/2 5n -/3n f m 2 (j 2 +k 2 )-Nm(j+k))/2 -jm -km km ^ 

0=0 j,k=0 

Despite the non-locality of the transformation, the Hamiltonian still remains local when expressed in 
terms of the parafermionic variables. The only complication arises for periodic boundary conditions 
since the naive application of the Jordan Wigner transformation generates a non-local term in this 
case. In other words, the correct Hamiltonian for the periodic parafermionic chain needs to be different 
from that of the periodic spin model, see Section 14.3.21 for details. 

In the following we will use a basis of parity eigenvectors in the space of ground states which is 
more natural to work with in the parafermionic picture: 

m— 1 

|0 a > := J>^"|^) , (4.3) 

/3=0 

so that 

<rf"|0 a ) = \<t> a+ p) , rt m \^a) = \<f>a) , P\M=u an \M- (4-4) 



4.2 Quantum phases 

In Section [3] we characterized SSB phases of Z^r spin chains. However, a moment of thought reveals 
that the local order parameter a an of the spin chain becomes non-local when expressed in terms 
of parafermionic variables. Since physical perturbations are assumed to be local, this observation 
suggests that there is no way to spontaneously break the symmetry in the parafermionic chain, and 
that the phases are characterized by topological order instead of ferromagnetic order. In this case, the 
states emerging at the macroscopic level would be those with definite Z^v parity. In the fermionic case 
N = 2, for example, it is well-established that there is no local order parameter, and correspondingly 
no local perturbation which can break the symmetry. 

While valid for a considerable subset of symmetry breaking patterns, the previous expectation 
fails to be true in general. This exactly happens in the phases labeled by a divisor n such that 
g = gcd(n,m) > 1. Parafermionic parity can be broken in these cases and, depending on the precise 
values of g and n, this leads to phases which, in addition to ferromagnetic order, might as well exhibit 
topological order. 

In order to substantiate our claim, let us briefly recall our results from Section 13.21 For g > 1 we 
have pointed out the existence of a second set of order parameters (|3.1ip . which when rewritten in 
terms of parafermions become local. More precisely, they correspond exactly to mutually commuting 
powers of the parafermions, 

„pN/g p'N/g 



(with p,p' e Z g ). (4.5) 

Note that, more generally, 7j JV/y commutes with 7"™ and 7^ m for every i,j, so that: 

\f /9 ,h t , l+1 ] =0. (4.6) 

In order to understand the physical implications of the order parameters let us then consider a 
Z^v — > /^g symmetry breaking perturbation of the original parafermionic Hamiltonian, 

H (->• H + hJ2^ N/9 + h-c.j (with peZg- {0}) . (4.7) 
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We suppose open boundary conditions, for which the Hamiltonians for the spins and the parafermions 
are identical. The extra term in the Hamiltonian induces the following action on the ground states, 

n-\ \pN/g 

t£3v«> = f n J °f /9 \^ = " apN/9 \^) > ( 4 - 8 ) 

and analogously for j^i, with an extra a-independent phase. As a consequence, the SSB induced by 
the perturbation (|4.7f) singles out a subset of the m ground states, exactly those labeled by an element 
a such that 

vN 

a — = (mod N) . (4.9) 
9 

It is evident that this equation has rh = m/g solutions 

a€{0,g,...,gm-l} = gZrh- (4.10) 

We believe that the degeneracy determined by the number of the remaining ground states cannot be 
lifted by SSB, and should therefore be regarded as topological. (If however m = g no such degeneracy 
is present, as e.g. when n = 2, N = 4, and the system is in a purely symmetry breaking phase.) 

Our previous interpretation receives additional support when considering the system with periodic 
boundary conditions. Namely, in that case the ground state degeneracy due to topologically protected 
edge modes is expected to disappear while the degeneracy due to symmetry breaking should persist. 
The degeneracy due to symmetry breaking is g, and therefore we expect a closed parafermionic chain 
to have a 5-fold degenerate ground state. Whenever g < m our arguments implies the existence of 
topologically protected edge modes. We illustrate this phenomenology in the case of Zs parafermions 
in figure [TJ Our expectations will be verified in Section T4.3I using the explicit realization of the phases. 
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n: l 2 4 8 n : 1 2 4 8 

Figure 1: (Color on-line) Degeneracies of the ground states as a function of the divisor n labeling the 
phases of Zs parafermions for both open (left) and closed (right) boundary conditions. 



4.3 Realization of the phases 

In this section we will provide an explicit realization of all parafermionic phases discussed in the 
previous section, thereby putting our considerations on a concrete basis. It turns out to be convenient 
to discuss open and closed boundary conditions separately. 

4.3.1 Open boundaries 

The Hamiltonian for the open chain is 

L-i 

i=l 
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where h^i+i is given in (|4.2p . H op has m degenerate ground states \<f> a ) ,a G Z m , see Eq. (|4.3p . 
When the phases are purely topological, namely when g := gcd(n,m) = 1, the degeneracy can be 
understood in terms of zero modes (operators commuting with Hamiltonian H op ) 7™ and 7^ which 
are localized at the boundary. Indeed, acting with powers of 7" (or j% L ) on a given ground state, say 
\4>o), produces the other m — 1 ground states distinguished by parity eigenvalues. In the generic case 
when g > 1, the ground state degeneracy is still due to zero modes, but now some of the m zero modes 
are no longer localized at the boundary. From Eq. (|4.6p we see that \^ N/9 ,H op ] = for every j, and 



acting with ^ N ^ 9 on 1 0o) will produce a state proportional to \<fi pm / g ), so that these zero modes are 
not located at the edge. These zero modes are precisely those responsible for the degeneracy due to 
symmetry breaking. In the particular case of purely topological phase labeled by n = 1 with m = N 
degenerate ground states, the question of robustness of the edge modes away from the finitely tuned 
Hamiltonians was raised and answered in [18]. We expect that all the zero modes responsible for 
topological degeneracy we found will also stay localized at the boundary in the whole phase. 

It has been pointed out in [U1[3B] that the quantum phases not distinguished by symmetry breaking 
(pure topological phases) correspond to different ways to "fractionalize" the symmetry operators on 
the space %i ow of low energy states of the open chain (or low entanglement energy states of the closed 
chain). According to these ideas, parity should be represented as the operator P = P\P T on 7^i ow 
in our case, with P\ and P r acting non-trivially only at the left and right boundaries, respectively. 
Denoting by x the Z^r degree of Pi (and by — x that of P r ), consistency of the factorization requires 
that [36] 



PP, 



P,P = P,P r P, 



l 



UJ 



pp. 



(4.12) 



In other words, the degree of Pi must satisfy the equation 

x(x - 1) = (mod N) . 



(4.13) 



The number of solutions to this equation with x € Z^r is given by the number of divisors n of N such 
that n and N/n are co-prime (see below). 

In what follows we will connect our previous considerations from Section 14.21 to the prediction 
obtained from symmetry fractionalization. Equation (j4.13j) selects out of the possible phases a distin- 
guished subset for which a low energy representation of parity exists. In those cases, the fractionalized 
symmetry operator can be constructed explicitly for our Hamiltonians (|4,lip . for which %i ow coin- 
cides with the ground state manifold. Since the edge modes are low energy operators localized at the 
boundary we make the ansatz 

(4-14) 



p = w r(iV+r)/2 7 r^-r 

The operator P satisfies the defining property of parity P N = 
space of ground states as 



1, is of degree zero, and it acts on the 



m— 1 



P\<j> a )=al(alyP r \<p a ) = u> anr J2 



0an a r\ 



0=0 



rn—l 



n-1 



n-1 



UJ 



£ ujP an u^ +mk > |/3 + mk) <g> 



+ mj) . 



(4.15) 
(4.16) 



0=0 



k=0 



3=0 



If we demand that P acts as P on the space of ground states, this requires that mkr = (mod N) 
for every k, that is r = ne, e € Z m . In this case, we have: 



P\K) 



UJ 



(4.17) 
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and e must be determined by the equation an(ne — 1) = (mod N) for every a € Z m . Clearly, the 
case a = 1 implies the others, so that we have to find an e such that 

n(ne-l) = (mod N) . (4.18) 

Alternatively, we can rewrite the previous equation as 

ne = 1 + pm with pGZ. (4.19) 

This equation has no solution if g := gcd (n,m) > 1. Indeed, in that case the left hand side of 
(ne — pm)/g = 1/g is an integer while the right hand side is not. On the other hand, if g = 1 there 
exists a unique solution!! This result is consistent with Section 14.21 where purely topological phases 
were predicted to occur only if g = 1. 

What is left to do is to show that we have exhausted all the purely topological phases predicted 
by Eq. (|4.13p . The operator P we have constructed, is obviously of the form P = PiP r with degree 
of P[ equal to ne, and Eq. (|4.13p is clearly implied by (|4.18p . We now prove that, furthermore, every 
solution of x(x — 1) = (mod N) defines a divisor n with gcd(n, N/n) = 1, so that solutions of 
Eq. (|4.13p and (|4.18p are in bijection. The proof is easily carried out. Setting n = gcd(x,N), simple 
properties of gcd give: 

? cd(n, N/n) = gcd ( g cd(x, N), ^J^) = ^gfe^^ (4.20) 

(4.21) 

(4.22) 



*cd(x,N) J gcd(x,N) 
gcd(gcd(x 2 , N 2 ),N) _ gcd(x 2 ,gcd(N 2 ,N)) 

gcd(x, N) gcd(x, N) 

gcd(x 2 (mod N),N) gcd(x,N) 



gcd(x,N) gcd(x,N) 

Therefore all the solutions of ()4.13j) can be expressed by those of (|4.18p in the form x = ne. 

We conclude with some more technical notes. Divisors n of N such that gcd(n, N/n) = 1 are called 
unitary. In the prime factorization iV = Hi??*, n = YIpT ^ s a unitary divisor if each c, is or Oj. For 
example the divisors of N = 4 are 1,2,4 but only 1, 4 are unitary. The number of parafermionic phases 
characterized by effective parity (exhibiting then pure topological order) is given by the number of 
unitary divisors. For N = 1, 2, 3, . . . their numbers are 

1, 2, 2, 2, 2, 4, 2, 2, 2, 4, 2, 4, 2, 4, 4, 2, 2, 4, 2, 4, 4, 4, 2, 4, 2, 4, 2, 4, 2, 8, . . . (Sloane's IA034444p . 

These numbers can also be identified with 2 q , where q is the number of different primes dividing N. 



4.3.2 Closed boundaries 

The parafermionic closed chain is defined by the local Hamiltonian 

Hpei = H op ~ h L l . (4.23) 

In specifying 1 we need to pay attention to the order of indices since parafermions on different sites 

do not commute. The substitution (2i,2i + 1) — > (2L,1) in the term 7<y 7^4.1 °f Eq. ([4.2p inverts 
the order of indices, and two different Hamiltonians are obtained depending on whether we commute 
the parafermions before or after the replacement of indices. Indeed performing the above substitution 
and then commuting the parafermions gives ^L^li^ 11 = UJ ^ n ^ li^ n l2L • ^ instead we commute the 

3 Uniqueness is easy to show since if e\ = (1 +pirn)/n and = (1 +P2m)/n are two solutions, then £1—62 = 
(pi — P2)m/n must be an integer, implying (pi — P2) is a multiple of n and so ei — €2 = (mod m). 
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parafermions in Eq. f|4.2[) before replacing the indices we will get uj~^ n ^ 2 ^^^l- The ambiguous 
extra phase produced by inverting the operations signals an ambiguity in the definition of h^i- The 
correct definition of the coupling instead keeps the order of labels and implements the replacement 
(2i,2i + 1) — > (1,2L) in J^^i+i- Following this prescription we get 



h P L1 = ^Y1 ^ n{N+Pn)/2 l P 1 n ^l n x Y, ^( m2 ^ 2+fe2 )^ Wm ^ +fc ))/^ 2 -^ l7 ^ 7 r fcm 7 2 fcm . (4.24) 

(3=0 j,k=0 

When rewritten in terms of spins, this operator turns out to be non-local. More precisely, it contains 
a twisting by the inverse global parity P': 

m— 1 n—1 

h it = ^ E (^wir e ^ m ■ ( 4 - 25 ) 

(3=0 j,k=0 

While completely equivalent when regarded with open boundary conditions, the spin model and the 
parafermionic chain thus exhibit essential differences when considered on a closed ring. 

In order to understand the physical implications of our previous comment, we now study the 
spectrum of the periodic chain, with particular attention to the ground state and its degeneracy. Recall 
that H op has m ground states which are given by \4> a ) as defined in Eq. (|4.3p . Since [H op , J = 0, 
the operators H op and h^ jl admit a common basis of eigenvectors, and our problem simply amounts 
to diagonalizing the projector hf il on the space spanned by the vectors \4> a )- 

A brief calculation shows that the operator l acts on the ground states of the open chain as 



^ m— 1 

- E = 5an =0 ( mod m)\(f>a) ■ (4.26) 



(3=0 



Put differently, the state \(f> a ) remains in the ground state manifold of the closed parafermionic chain 
if and only if 

an = (mod m) . (4.27) 

In order to enumerate the solutions to this equation, let us first introduce the symbols n = n/g and 
m = m/g where we used the previous abbreviation g = gcd(n, m). With this notation, we can write 
the solutions of Eq. ([4.27|) as 

m 

a = p— with p G nZ . (4.28) 
n 

Hence, the possible values of p are p = 0, n, . . . , h{g — 1), so that there are g distinct solutions. 

It may easily be shown that the previous g states exhaust the ground space manifold of the periodic 
Hamiltonian (|4.23p . Note that we can associate the degeneracy of the closed chain to parafermionic 

zero modes as in the case of the open chain. Indeed Eq. (|4.6p implies that 7? N ^ 9 , = 0, and 

acting with powers of on a reference ground state will produce g mutually orthogonal ground 
states with different parity eigenvalues. 

Our analysis implies that the closed parafermionic chains have a unique ground state for n being a 
unitary divisor, while there exists a g-fold degenerate ground state when n fails to be a unitary divisor. 
In the former case, there is no spontaneous symmetry breaking and the ground state degeneracy in the 
open chain may be completely attributed to the non-trivial topology of the system. On the other hand, 
in the latter case, the ground state degeneracy is only partially lifted when making the transition from 
open to closed boundary conditions. In that case, we observe a combination of non-trivial topology 
and spontaneous symmetry breaking if m/g > 1. For g = m there is no non-trivial topology and only 
symmetry breaking. These results are in complete agreement with our expectations from Section 14.21 
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5 Conclusions and outlook 



In the current paper we have classified the massive phases of Zjy parafermionic chains. This was 
achieved by constructing all gapped symmetry broken phases in the dual quantum spin chains. The 
number of distinct phases and their specific characteristics depend crucially on the number-theoretic 
properties of N. First of all, it is easy to construct one SSB phase of the spin model for each divisor 
n of N. These phases are characterized by the (co-) existence of local and non-local order parameters. 
Since Zjy does not lead to symmetry protected topological phases [12 [29], the enumeration of divisors 
exhausts all possibilities for gapped phases. 

In a second step we have then interpreted the resulting ground states and Hamiltonians from a 
parafermionic perspective. Of course, this leads to the same number of phases. However, in contrast 
to the spin model the parafermionic phases can now exhibit features of both, symmetry breaking and 
non-trivial topology. Whether the resulting parafermionic phase should be interpreted as topological 
first of all depends on whether n and N/n still have common divisors. A parafermionic phase is purely 
topological if n is a unitary divisor of N. By definition, this is the case if n and the quotient N/n are 
co-prime or, in our notation, if g = gcd(n, N/n) = 1. The number of unitary divisors is known to be 2 q 
where q is the number of distinct prime factors in N. In contrast, a parafermionic phase described by 
n < N not being a unitary divisor (i.e. g > 1) definitively exhibits spontaneous symmetry breaking. 
If, in addition, gn < N then the phase is also protected by topology. Topology protected phases in 
open systems exhibit gapless parafermionic edge modes. 

While our current results are only concerned with the intrinsic "Ln symmetry of parafermionic 
chains, it would be interesting to extend our analysis to systems which are required to respect ad- 
ditional symmetries such as inversion or time-reversal symmetry. Alternatively, one could also add 
internal degrees of freedom which transform under a discrete or continuous symmetry group. In all 
these cases one is lead to so-called symmetry protected topological phases which, for bosonic systems, 
have been fully classified in Refs. [311271 [29], It should be noted that models with continuous symmetry 
lead to a few peculiarities that have been addressed in Ref. [37] . 

For fermionic chains it is known that the inclusion of additional symmetries has profound conse- 
quences. Most importantly, the restriction to time-reversal invariant systems enhances the group of 
topological invariants to Z or Zg, respectively, depending on whether interactions are allowed |381126[ [6] 
or not [U [5] . A general mathematical formalism for treating fermionic symmetry protected topolog- 
ical systems has been developed in Ref. [39] and it appears feasible to extend it to the Zjy-graded 
generalization of the Clifford algebra |28j needed to describe parafermions. 

Finally, the classification scheme used in this paper is, of course, not limited to Z^v clock models. 
It should rather be applicable to any model of statistical physics which exhibits a Kramers- Wannier 
type of duality. While abelian dualities are well-explored (see, e.g., [30]) and close to our exposition, 
recent progress on non-abelian dualities [41J might offer a vast playground for the exploration of new 
topological phases. 

Needless to say, knowledge about the classification of topological phases, symmetry protected 
or not, and a representative Hamiltonian for each phase should merely be the starting point for a 
more thorough investigation. Indeed, the ultimate goal should be to gain control over the complete 
phases diagram of physically realistic Hamiltonians, including the location and the nature of phase 
transitions. For the Majorana chain with interactions such an analysis was provided in [22], based 
on the known phase diagram of the dual anisotropic next-nearest-neighbor Ising (ANNNI) model. It 
would be interesting to extend this analysis to the case of Z^ spin models. Also, with regard to the 
potential physical realization and detection of parafermions one will need to understand the effects of 
adding disorder to the couplings. 

* * * Note * * * 

While in the process of completing this work, the paper [36J appeared which has considerable overlap 
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with our own results. We gratefully acknowledge useful discussions with J. Motruk and A. Turner 
about aspects of entanglement and symmetry fractionalization in parafermionic chains prior to the 
publication of their paper. 
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